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We present a theory for Raman scattering on 2D quantum antiferromagnets. The microscopic 
Fleury-Loudon Hamiltonian is expressed in terms of an effective 0(3) - model. Well within the Neel 
ordered phase, the Raman spectrum contains a two-magnon and a two-Higgs contribution, which 
are calculated diagramatically. The vertex functions for both the Higgs and magnon contributions 
are determined from a numerical solution of the corresponding Bethe-Salpeter equation. Due to the 
momentum dependence of the Raman vertex in the relevant B\ g + Eig symmetry, the contribution 
from the Higgs mode is strongly suppressed. Except for intermediate values of the Higgs mass, it does 
not show up as separate peak in the spectrum but gives rise to a broad continuum above the dominant 
contribution from two-magnon excitations. The latter give rise to a broad, asymmetric peak at 
to ~ 2.44 J, which is a result of magnon-magnon interactions mediated by the Higgs mode. The 
full Raman spectrum is determined completely by the antiferromagnetic exchange coupling J and a 
dimensionless Higgs mass. Experimental Raman spectra of undoped cuprates turn out to be in very 
good agreement with the theory only with inclusion of the Higgs contribution. They thus provide a 
clear signature of the presence of a Higgs mode in spin one-half 2D quantum antiferromagnets. 


I. INTRODUCTION 


The recent observation of a scalar Higgs particle as 
the final missing ingredient of the standard model in 
particle physics 1,2 has led to a resurgence of interest in 
the question whether a well defined excitation associated 
with an amplitude mode is present and observable in 
condensed matter systems with a broken continuous 
symmetry 3 . A case in question are neutral superfluids, 
in which a global U( 1 ) symmetry is broken as the 
simplest example of continuous symmetry breaking. The 
dynamics of the superfluid order parameter, however, 
is generically of the Gross-Pitaevskii type, i.e. first 
order in time. Despite the Mexican hat structure of 
the effective potential, there is thus only a Goldstone 
but no Higgs mode. Indeed, a necessary condition for a 
Higgs mode in the non-relativistic context of condensed 
matter systems is an emergent Lorentz invariance. This 
appears, for instance, in the superfluid phase of ultracold 
Bosons in an optical lattice at integer filling, where the 
particle-hole symmetry close to the transition to a Mott 
-insulating phase gives rise to an effectively Lorentz 
invariant theory 4,5 . In this regime, a Higgs mode has 
indeed been observed from the absorbed energy in 
response to shaking the optical lattice 6 . The associated 
effective Higgs mass vanishes in a continuous manner 
at the transition to a Mott insulator, essentially like a 
mirror image of the gapped particle-hole excitations in 
the incompressible phase. A Higgs mode has also been 
observed via neutron scattering in the antiferromag- 
netically ordered phase of TICUCI3 at high pressure 7,8 . 
In both cases, the evolution of the Higgs mode can be 
followed into the regime where the order parameter 
vanishes near a quantum critical point. Much earlier, 
in superconductors, the presence of an amplitude mode 
has been inferred by Sooryakumar and Klein from the 
observation of a distinct peak in Raman scattering from 


NbSe2 in a regime, where superconductivity coexists 
with a charge density wave 9 . In this special situation, 
the superconducting order parameter directly couples to 
light because the modulation of the charge density wave 
by the electromagnetic field in the A\ g symmetry of the 
Raman response leads to a periodic modulation of the 
density of states at the Fermi surface and thus of the 
superconducting gap 10 . The amplitude mode associated 
with oscillations of the gap magnitude can thus be 
directly seen in Raman scattering. This interpretation 
of the experiment has recently been confirmed within a 
detailed theoretical analysis by Cea and Benfatto 11 and 
also in an independent experiment by Measson et.al. 12 . 

In our present paper, we analyze the possibility to 
observe a Higgs mode in quantum spin models whose 
ground state exhibits antiferromagnetic order in two di¬ 
mensions via Raman scattering. This problem is of in¬ 
terest for at least two different reasons: First of all, 
there are experimental data on clean samples of undoped 
cuprates 13,14 for which no theoretical model and under¬ 
standing seems to have been given so far. Second, due 
to its relevance in the parent compounds of high tem¬ 
perature superconductors, quantum antiferromagnets in 
two dimensions are among the most intensively studied 
examples of a broken continuous symmetry in condensed 
matter. In particular, the issue of a well defined Higgs 
excitation is nontrivial in two dimensions ( 2 D). Indeed, 
in this case the phase space for the decay of a Higgs mode 
into two Goldstone modes is very large. As a result, the 
Higgs spectrum at zero wave vector extends down to zero 
energy, diverging with a universal power law ~ 1/w 15,16 
which has in fact been observed in neutron scattering 
from the Neel state of undoped La2Cu04 4 '. To leading 
order, this low energy singularity is due to two-magnon 
processes. The presence of a large weight in the spec¬ 
trum of the amplitude mode at low energies suggests that 
the Higgs mode is overdamped in 2 D. As emphasized by 
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Podolsky, Auerbach and Arovas 18 , this is not the case in 
general, however. In particular, response functions which 
- unlike neutron scattering - couple to the square of the 
order parameter are expected to show a clear Higgs peak 
even in two dimensions because for them the low fre¬ 
quency response is suppressed by a factor w 3 . Precisely 
this type of response is measured in the lattice shaking 
experiment near the superfluid to Mott insulator tran¬ 
sition performed by Endres etal 6 . As a result, a sharp 
peak due to a Higgs excitation appears even in 2D, as 
supported by a detailed theoretical analysis of this prob¬ 
lem, both numerically 5,19 ’ 20 and analytically 21 . In par¬ 
ticular, the Higgs mode remains well defined even near 
the quantum critical point because both its energy and 
width vanish with the same power in the distance from 
the critical point 5,19 ” 22 . 


II. RAMAN SPECTRUM 


A. Effective light scattering operator 


The differential cross section for Raman scattering is 
given by 26 


cPcr 

d£LdhJ 


= hr 2 0 ^R(i^ f ) 


(1) 


where ro = e 2 /mc 2 is the Thompson radius and R(i —> /) 
is the transition rate from the initial photon state |k;, e,) 
to the final state |k/,ey). The rate is determined by 
Fermi’s golden rule 25 ’ 2 ' 


R(i^f) = ^Y. e ~ PEl \ M F,i\ 2 S( E F- E i-hco) (2) 

I,F 


Regarding the Neel ordered state of 2D quantum an- 
tiferromagnets, Raman scattering seems to be ideally 
suited to observe the amplitude mode since it couples 
to the square of the order parameter 18 . Experimentally, 
however, no indications for a Higgs mode are found 13 . In¬ 
stead, there is a broad, asymmetric peak which appears 
to be due to a strongly renormalized two-magnon excita¬ 
tion. As will be shown below, this observation can be ex¬ 
plained both in qualitative and quantitative terms within 
a rather simple 0(3) - model of the Neel ordered phase. 
In particular, the absence of a separate peak associated 
with the Higgs mode is a consequence of the fact that 
the leading order contributions from the amplitude mode 
vanish because of the nontrivial Raman vertex in the rel¬ 
evant Big symmetry. Nevertheless, the Higgs mode turns 
out to have a crucial effect on the detailed form of the Ra¬ 
man spectrum. In the dominant two-magnon response, 
it shows up by mediating the magnon-magnon interac¬ 
tions, which lead both to a downward shift of the peak 
compared to a simple spin-wave analysis 23,24 and a char¬ 
acteristic asymmetric line shape. Moreover, the Higgs 
mode gives rise to a broad continuum above the two- 
magnon peak. Both features are found in experimental 
spectra for La 2 Cu 04 14 , YBa 2 Cu 3 C> 6 13 and also in - still 
unpublished - data on Sr 2 Ir 04 25 . Using the well known 
values for the antiferromagnetic exchange coupling J in 
these systems, the data turn out to be in very good agree¬ 
ment with the theory up to frequencies of order 6 J and 
beyond. As will be shown below, the agreement between 
experiment and theory relies crucially on whether the 
Higgs mode contribution to the Raman response is in¬ 
cluded or not. Raman scattering can therefore indeed 
be used to probe the existence of a Higgs mode in 2D 
quantum antiferromagnets. 


where 1 1), |.F) are the initial, final state of the sample, M 
is the effective light scattering operator and to = uii — uj / 
is the frequency transfer. For the specific example of 
the Mott-insulating, antiferromagnetic phase of undoped 
cuprates, an appropriate microscopic Hamiltonian is the 
nearest neighbor Heisenberg model on a square lattice 
which arises from a superexchange mechanism 28 . The as¬ 
sociated effective light scattering operator is then given 
by the Fleury-Loudon Hamiltonian in leading order of a 
moment expansion 2 ' ,29,3 ° Performing the common sym¬ 
metry decomposition of the Raman scattering opera¬ 
tor 26 , it turns out that only the Hi g -mode is Raman ac¬ 
tive because the operators of the other symmetry modes 
commute with the Heisenberg Hamiltonian. Restricting 
the interaction to this mode, the reduced effective Hamil¬ 
tonian is given by 

H fl = 2BP(e i ,e / ) ^ (s Xj ■ S x .+* - S x . ■ S x .+y) 

3 

(3) 

with P(ei,ef) = efe^ — e\e v ^ and B ~ t 2 /(U — hcoi). 
As has been shown by Shastry and Shraiman 2 ', there 
is also a response in other symmetry modes associated 
with ring exchanges on a plaquette or the fluctuations 
of a chiral spin operator 31 . They require an extension 
of the microscopic Hamiltonian beyond the Heisenberg 
model with nearest neighbor exchange only and thus lead 
to corrections to the simple form of the Fleury-Loudon 
Hamiltonian (3) . In our present work, we do not consider 
such extensions. In fact, as will be shown below, already 
the leading order Hamiltonian (3), which applies to 2D 
quantum antiferromagnets with a well defined Neel order, 
properly accounts for the essential features of the Raman 
spectrum. Apart from microscopic prefactors, the tran¬ 
sition rate in Eq. (2) can then be written as the Fourier 
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transform of the van Hove function S(t) of L/fl 


f) = ^9(ki)g(k f ) J dt e iut (H FF (t)HM0)) T 

—oo 

= ^2 9 (ki)g(kf)S(u}) (4) 

The resulting Raman response function S(ui) is con¬ 
nected to an associated spectral function \R am an( w ) by 
the fluctuation-dissipation theorem. At zero tempera¬ 
ture, this takes the simple form 


SM = 2fiXRaman(w)0(w) = 2/l ImX'R a man(w + i0)0(u}) 

(5) 

In practice, the spectral function will be determined from 
the imaginary time ordered correlation function 


XRa 


i{t) = h (Tt 




( 6 ) 


by analytic continuation to real frequencies. 


B. Field theoretic formalism 

In order to disentangle the contributions from magnons 
and the Higgs mode to the Raman spectrum, it is both 
necessary and convenient to replace the spin-operators 
of the microscopic Heisenberg model by a coarse grained 
description in which the relevant excitations appear more 
directly. Such an effective description is provided by the 
0(3) nonlinear cr-model (NLcrM), which is the continuum 
field theory of the antiferromagnetic Heisenberg model as 
used by Chakravarty, Halperin and Nelson 32 . The asso¬ 
ciated effective action can be written in an apparently 
relativistic invariant form 



with a three-vector (x* 1 ) = (c s r, x) and spin wave veloc¬ 
ity c s . Apart from c s and the momentum cutoff A, the 
dimensionless coupling constant g = gA of the nonlinear 
cr-model contains the spin stiffness p s , which is of the or¬ 
der of the exchange coupling constant of the underlying 
Heisenberg model. The slowly varying order parameter 
field n(a;) is a three component field which describes de¬ 
viations from perfect Neel order. It is normalized by 
n 2 = 1. On the microscopic level, the physical meaning 
of the field n becomes evident by the mapping 

fel (_i) s g n ( x 3) n ( Xj .)y/f _ a 4 L 2( Xj .) + a 2 L(x,) 

( 8 ) 

between spin operators in the semiclassical limit S 3> 1 
and bosonic fields originally due to Haldane 33,34 . Here a 
is the lattice constant and the angular momentum field 
L is constrained by n-L = 0, a 4 L 2 -C 1. Using Eq. (8) to 


express the spin operators in the Fleury-Loudon Hamil¬ 
tonian (3) in terms of the order parameter field, it turns 
out that the Raman spectrum requires to determine the 
connected correlation function of the following operator 
in the NLcrM 

O(t) = 2 B H 2 S 2 P(e.i,ef) J d 2 x cr|)<9jn( X , r) • 9jn(x,r) 

= 4 BH 2 S 2 P{e l7 e f ) J ^^ 7 (k) |n(k, r)| 2 (9) 

Here 7 (k) = \ {k 2 — k 2 ) is the continuum form of the Big- 
symmetry factor. Unfortunately, the combination of non¬ 
trivial vertices in the NLcrM which appear beyond a sim¬ 
ple Gaussian approximation due to the constraint n 2 = 1 
and the complicated form of the operator (9) make it very 
hard to calculate the Raman spectrum directly within the 
NLcrM. Moreover, it is rather difficult to disentangle the 
contributions which arise from magnons or the Higgs ex¬ 
citation in the standard decomposition n = (tv, \/l — 7r 2 ) 
which only involves the two-component field tv associated 
with the magnons. It is therefore more convenient to use 
a soft-spin description of the order parameter within a 
linear 0(3) cr-model, whose effective action 


sm 




( 10 ) 


contains the bare mass mo of the amplitude mode as an 
additional parameter. It determines the stiffness for lon¬ 
gitudinal fluctuations in the Neel ordered phase which 
arise from the fact that for quantum spins the operator 
Sqaf which determines the long range antiferromagnetic 
order at the wave vector Qaf has finite fluctuations along 
the ^-direction 35 In general, the coupling constant g and 
to o are independent parameters. It is only in the scaling 
regime close to the quantum critical point g = g c for the 
loss of Neel order, where the value of the bare Higgs mass 
is irrelevant and where the 0(3) cr-model is completely 
equivalent to the NLcrM 21 . In fact, as will be shown be¬ 
low, the dominant features of the Raman spectrum are 
not very sensitive to the value of the Higgs mass Too even 
in the opposite regime g <C g c of well defined Neel order. 
They are determined essentially by the spin stiffness or 
the equivalent Heisenberg exchange energy J, which is 
known rather accurately from neutron scattering data. 
An intuitive argument for why it is possible to replace the 
fixed length NLcrM with a soft-spin model relies on view¬ 
ing the linear cr-model as an effective low-energy theory 
of the NLcrM, where short wave-length fluctuations have 
been integrated out 36 . The Neel-field n(x) is therefore 
averaged over some domain in real space and the result¬ 
ing averaged field $ is no longer constrained to have a 
fixed length. Instead, it is subject to a potential U(3>), 
which exhibits a minimum at a finite vacuum expectation 
value (VEV), which has been chosen to be (|3?|) = x/3 in 
Eq. (10). The magnitude of the fluctuations around the 
minimum is controlled by the dimensionless mass param¬ 
eter mo = toq/A, which is expected to be large compared 
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to one in the semiclassical limit S 1, where fluctua¬ 
tions of the order parameter magnitude vanish. The dy¬ 
namics of the averaged field is fixed by the condition that 
the magnon dispersion Wfc = c s |k| is linear in momentum, 
giving rise to an effectively Lorentz invariant action. 

In order to calculate the Raman spectrum within the 
linear cx-model (10), we use the standard parametrization 
for the field $ in the symmetry broken phase: 

$ = ^7r, r \/3 + a'j . (11) 

The parameter r = r(g, A) incorporates the renormaliza¬ 
tion of the VEV by quantum fluctuations and has to be 
determined from the condition (cr) = 0. The components 
of 7r are the two massless antiferromagnetic magnons, 
while a denotes the massive amplitude or Higgs mode. 
Expressed in terms of the relativistic three-momentum 
k = (oj/c s , k), the free propagators for the fields are 


which involve a two-magnon susceptibility xm and a 
two-Higgs susceptibility Xh- They are determined by 
the correlation functions of two Magnon and two Higgs 
operators Om and Oh , respectively. In addition, an in¬ 
terference susceptibility xint arises which consists of the 
two mixed correlation functions. A similar structure has 
been found by Canali and Girvin 38 , who have calculated 
the Raman spectrum of undoped cuprates by means of a 
Dyson-Maleev transformation of the underlying antifer¬ 
romagnetic Heisenberg model. Compared to our present 
formulation, this approach does not allow to separate 
the magnon and Higgs contributions to the spectrum, 
with the latter appearing as part of a four-magnon 
term. Moreover, the magnon-magnon interaction is 
replaced by an instantaneous one, thus missing the 
retardation effects associated with the Higgs-mediated 
interaction discussed below which is treated properly 
by our numerical solution of the Bethe-Salpeter equation. 


G°^(fc) = ('Ki(k) / K j (—k)) = 5ijG® n (k) = Sij-jL (12) 
G° CT (fc) = Hk)a(-k)) = (13) 

/C TRq 


These propagators will change with increasing strength 
g of the quantum fluctations. In particular, the exact 
propagator 


G arJ (k) = 


TOo 


k 2 


0 k 2 


]f 24|fc| 


(14) 


To leading order in the coupling g , only the two- 
magnon and two-Higgs susceptibilities are nonzero while 
the interference susceptibility contributes only at order 
g 3 . The corresponding Feynman diagrams in FIG. 1 
are just a two-magnon, respectively a two-Higgs bubble 
dressed with two Raman symmetry factors q(k). Drop¬ 
ping the polarization factor P(e i7 ey), the two-magnon 
response to lowest order in g is given by 

n 2 R 2 fi 3 

XmM = -—Xo -w 3 6»(2c s A — u) + 0(g 3 ) . (19) 

12 c s 


of the Higgs mode acquires a contribution ~ \/\k\ due 
to the decay into Goldstone modes. Instead of a simple 
pole, the Higgs propogator therefore exhibits a branch 
cut which starts at the threshold uif- = c s |k| for the ex¬ 
citation of magnons 16 . This low energy singularity first 
appears at order g 2 and is expected to remain valid to 
all orders, i.e. G aa (k —► 0) is proportional to l/\k\ at 
any finite value of g. Similarly, by the Goldstone the¬ 
orem, the magnons are always massless. The structure 
G- 7I T T (k) = Zg/k 2 therefore remains exact with a finite, 
renormalized coupling constant Zg throughout the Neel 
ordered phase 37 . Inserting Eq. (11) into Eq. (9) we ob¬ 
tain 


O(t) = OmW+OhW 


O m (t) 


{d x nf 


0 H (r) oc / d 2 x ( d x a) 


(d v n) 2 

(d y a) 2 


(15) 

(16) 

(17) 


The operator to which light couples in the Raman re¬ 
sponse therefore involves the square of gradients of the 
magnon and Higgs field. As a result, the response al¬ 
ways involves at least two magnon or Higgs excitations. 
Taken together, the full Raman susceptibility is the sum 
of three different contributions 


XRamanM = Xm( W ) + XhM + XlntM (18) 


The spectrum is a pure power law ~ u> 3 with a sharp cut¬ 
off at twice the zone-boundary energy of a magnon. If we 
assume, that the Heisenberg model underlying the con¬ 
tinuum theory contains only a nearest-neighbor coupling 
J, this corresponds to 2c s A = 2nhJ. 

Regarding the contribution from longitudinal fluctu¬ 
ations of the order parameter, the spectrum to leading 
order in g is determined by the diagram in FIG. 1 b, 
which gives 


XhM = 


g 2 B 2 H 3 S 4 
12c, w 


(!) 


1 2 


- - c,m n 


6(uj — 2c s mo) 


9 (^2c s \JA 2 + TOq — + 0(g 3 ). 


( 20 ) 


It shows the expected threshold at twice the Higgs mass 
and peaks at the maximum energy of two Higgs excita¬ 
tions. Note that in contrast to neutron scattering, which 
couples linearly to the order parameter, the Raman spec¬ 
trum does not contain the longitudinal propagator G aa 
directly. As a result, no sharp peak due to the Higgs 
mode appears even at leading order. In the following 
section, we will show that the sharp features in both the 
two-magnon and the two-Higgs response, which depend 
quite sensitively on the precise form of the momentum 
cutoff and the specific value of the Higgs mass mo, are 
completely eliminated in a calculation which sums up the 
next-to-leading order diagrams in a consistent fashion. 
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FIG. 1. Leading order diagram for the two-magnon (a) and 
the two-Higgs susceptibility (b). The dashed line propagators 
are those of the bare magnons in Eq. (12), the full lines denote 
the Higgs propagator in Eq. (13). 


(a) 'X/X/X* ^ «W\/ 


(b) -wx* •—.—i «w\. 


(c) 'x/x/x* 


)<x/x/v 


FIG. 2. Next-to-leading order diagrams apart from self-energy 
insertions for the 2-magnon susceptibility \m 


C. Beyond leading order and Bethe-Salpeter 
equation 

The next-to-leading order diagrams, apart from self¬ 
energy insertions, for the two-magnon, two-Higgs and 
interference susceptibility are shown in figures FIG. 2, 
FIG. 3 and FIG. 4 respectively. To see whether these 
diagrams may lead to sharp features associated with the 
Higgs mode, we start with a qualitative discussion based 
on symmetries alone. One immediately notices, that the 
loop integrals in FIG. 2(a) decouple at vanishing external 
wave vector, where q = (to, 0). These diagrams therefore 
give no contribution to the Raman response because of 
the antisymmetry of the f?i s -symmetry factor y(k) under 
the exchange k x k y . 

Especially the vanishing of diagram FIG. 2(b) is im¬ 
portant to notice, as it would lead to a peak at the Higgs 
mass due to the pole in the Higgs propagator at least in 
leading order in g. The only non-vanishing diagram for 
the two-magnon susceptibility is FIG. 2(c), in which the 
Higgs line is independent of the frequency transfer. As a 
result, no sharp peak at the Higgs mass arises. For the 
two-Higgs and the interference susceptibilities the situa¬ 
tion is quite similar. Indeed, the diagrams (a) and (b) 
in FIG. 3 and FIG. 4 vanish due to the antisymmetry of 
the I?i g -symmetry factor under the exchange k x O k y . 
Concerning the interference contribution, it turns out, 
that up to order g 3 it gives a negligible contribution to 
the full Raman response. Thus, in the following, we only 


FIG. 3. Next-to-leading order diagrams apart from self-energy 
insertions for the 2-Higgs susceptibility xh 




FIG. 4. Next-to-leading order diagrams for the interference 
susceptibility xint- 


consider the two-magnon and two-Higgs contributions. 

In order to determine the Raman spectrum in quan¬ 
titative terms, we now derive a formally exact equa¬ 
tion for the two-magnon and the two-Higgs response. 
Based on the diagrammatic expansion shown in FIG. 5, 
the two-magnon contribution \m can be written in 
terms of the exact two-magnon Raman vertex function 
T (k, k', 12,12', w) in the following form 

/ rl 3 k I" rl 3 k' 

(2^ J ^( k )7(k')G..(k,I2 + W ) 

X G 7r7r (k, f2)r w (k, k', 12,12', w) . 

( 21 ) 

The vertex function obeys a Bethe-Salpeter equation 

r fj3 ur 

r w (k,k',I 2 ,I 2 ',a;) = 2(27r ) 3 6(k-k') + J ^V*(fc, k", u) 

x G nn ( k", 12" + w)G W7r (k", 12")r T (k", k', 12", fi» 

( 22 ) 

which is depicted diagrammatically in FIG. 6. To deter¬ 
mine the Raman spectrum, we do not need the full vertex 
function. Instead, it is sufficient to know the reduced one, 
which is defined as 

r fj3u 

r T (k,o, W ) = J ™ 7 (k , )r T (k,k , ,n,n , , w ) ( 23 ) 

Using this definition, the equations for the two-magnon 
susceptibility, Eq.’s (21) and (22), can be recast in the 
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FIG. 5. Diagrammatic series for the 2-magnon susceptibility. 
Magnon lines have to be interpreted as full propagators. The 
shaded area represents the full 2-magnon Raman vertex 


k. + u) k, 0, 4- u) k", Q" + u 



k, n k", Q" 


g —> Zg. Regarding the Higgs-propagator, the leading or¬ 
der result contains a simple pole at the bare Higgs mass 
mo but misses the l/\k |-singularity due to the decay of 
the Higgs into two Goldstone modes, see Eq. (14). The 
final form of the Bethe-Salpeter equation in ladder ap¬ 
proximation is given by 


r T Ladder (k,H, w )= 27 (k) + ^ J 0^GUk-k") 

X G° wn ( k", f]")G° w (k", f r + w )r“ r (k", n", w) . 

(28) 

This is still too complicated to be solved analytically, 
however it is amenable to an efficient numerical solution. 
Before, doing so, however, a similar equation will be 
derived for the two-Higgs susceptibility. 


FIG. 6. Diagrammatic series for the full 2-magnon Raman 
vertex, giving the Bethe-Salpeter equation (22). 

form 

/ d 3 k 

—^7(k)G W7r (k, n + w) 

xG„(k,!l)r,(k,S], w ) (24) 

r tj3ur 

r T (k,n,w) = 2y(k) + J j^v v (k,k J, ,u)G„ v ( k",n") 

x G™(k",B"+u;)r w (k",B",w) (25) 

Up to this point, there are no approximations involved. 
A solution of Eq. (25), however requires knowledge of 
both the exact magnon propagator 


The way to derive the integral equations for the two- 
Higgs susceptibility is completely analogous to the case of 
the two-Magnon susceptibility, one just needs to replace 
the magnon propagators by Higgs propagators. Hence we 
will just write down the resulting equations, containing 
already the reduced two-Higgs Raman vertex function 

^(k.n.w): 

/ f fi i. 

^p7(k)G CTCT (k, Q + u})G aa { k, 0) 

xU(k,Sl,w) (29) 

r /)3 u" 

r CT (k,H, W )= 27 (k)+ J ^V a (k,k J, ,u)G aa (k",tf') 
x G a<J (k",n" + co)r a (k'\n",oj) (30) 


Gnir{p) 


9 

p 2 - gTj^(p) 


(26) 


including all self-energy corrections via the exact magnon 
self-energy and the full irreducible two-magnon in¬ 
teraction V^(k, k", oj). 

For an explicit and quantitative solution, we truncate 
the irreducible magnon interaction at lowest order in g 1 
which amounts to a ladder approximation. Properly tak¬ 
ing into account combinatorial factors we find 

V ir (k,k",u) = 2-4.i(^§-) G® a (k — k”) (27) 

Thus, to lowest order, the interaction of magnons which is 
relevant for the B\ g symmetry is mediated by the Higgs 
mode. Inserting Eq. (27) into Eq. (25), we obtain the 
ladder Bethe-Salpeter equation, which corresponds to a 
consistent next-to-leading order approximation in the ir¬ 
reducible magnon-magnon interaction. For the magnon 
propagator, the zeroth order result G° n captures the 
exact low energy behavior ~ g/\k\ 2 due to the Gold- 
stone theorem apart from a finite renormalization factor 


The leading diagrams for the full irreducible two-Higgs 
interaction V a (k,k",u>) are depicted in FIG. 7. Truncat¬ 
ing the irreducible two-Higgs interaction again at leading 
order in g , we obtain 


V a (k, fc",w) = 2 • (2 • 3) 2 • l - G° aa (k ~ k") 

(31) 

Note that this has precisely the same form as the ir¬ 
reducible two-magnon interaction, differing only in the 
combinatorial and coupling prefactors. Inserting eq. (31) 
into Eq. (30) we obtain the Bethe-Salpeter equation for 
the reduced two-Higgs Raman vertex function 

r ff “ r (k,H, w ) = 27 (k)+ 3 p° J ^G° ff (fc-fc") 

X g° ct (k", B")G° CT (k", n" + w)r“ r (k",n",w) 

(32) 


which again needs to be solved numerically. 
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k, fl + lo k", 0" + u 


V CT (fc, k", co)= 


fc-fc" + 0(g 2 


k, n k", n" 


FIG. 7. Diagrammatic representation of the irreducible 2- 
Higgs interaction vertex. 


D. Numerical solution of the Bethe-Salpeter 
equation and Raman spectra 


For the numerical solution of the two independent 
Bethe-Salpeter equations (28) and (32), it is conve¬ 
nient to first write them in dimensionless form. This is 
achieved by introducing k = k/ A, mo = mo/A, g = gA, 
x(w) = x(cb)/A and f(k,k',w) = A _2 F(k, k', (b).The di¬ 
mensionless frequency Co = u>/(c s A) is measured in units 
of c s A = 7 tHJ ~ S. Note that the dimensionless cou¬ 
pling g = gA is related to the spin quantum number S 
by g = 27r/5, which is far from small in the relevant 
case 5=1/2 discussed below. In practice, the relation 
g = 27r/5, which arises from a semi-classical expansion 
valid for S> 1, does not enter our results for the Raman 
spectra in quantitative terms, except for determining the 
relative weight and frequency scales for different 5 in 
Figs. 8 and 9. For a given value of 5, our expansion to 
leading order in g is therefore valid provided g -C g c , he. 
the antiferromagnet is far away from a complete destruc¬ 
tion of Neel order by quantum fluctuations 32 . 

The explicit numerical solution of the Bethe-Salpeter 
equation employs the standard approach to integral equa¬ 
tions, namely discretizing the integrals and solving the 
resulting system of linear equations. For the dimension¬ 
less Brillouin-zone [— 1, l] 2 , we use 20 x 20 points. Since 
an increase up to a grid size of 30 x 30 gave only very 
small corrections, we can assume that the results have 
converged to the continuum limit. While the frequency 
integration has no intrinsic microscopic cut-off, it is re¬ 
stricted in practice to values below 4c s A. Contributions 
beyond that do not affect our results. We have varied 
the number of sampling points for the frequency integral 
between 20 and 40. Similar to the situation for the mo¬ 
mentum integration, the grid size had no influence on the 
results. 

In FIG. 8, the two-magnon susceptibility Xm(w) ob¬ 
tained from the numerical solution of the Bethe-Salpeter 
equation is shown for three different values of the dimen¬ 
sionless coupling g , corresponding to 5 = 1/2,1,2 and 
for an intermediate value of the Higgs mass in o = 0.5. 
To obtain the associated spectral function, which deter¬ 
mines the Raman spectrum in real frequencies, an ana¬ 
lytic continuation is required. Since our numerical data 
contain no noise, we use an Ansatz to fit the data which 
can then be continued in analytical form. Specifically, as 



FIG. 8. (Color online) Numerical results for the 2-magnon 
susceptibility as a function of imaginary frequency for S = 
1/2, S = 1 and 5 = 2. The solid curves are fits of Eq. (33). 


a fit model, we use 


X(w) 


«i + a 2 q 2 

1 + a 3 q 2 + a 4 g 3 + a 5 q 4 



(33) 


This Ansatz is motivated by our leading order results for 
the two-magnon susceptibility, which obeys xm ~ w~ 2 
for large Co and Xm ~ w 3 f° r small Co. The fitted curves 
are also shown in FIG. 8. Apparently they reproduce the 
numerical results very well with a variance which is only 
1.4 x 10" 6 . 

Based on the explicit form of the spectral function 
Eq. (33) in terms of the imaginary frequency variable 
q, the analytic continuation q —> — iCo is easily done. 
The resulting two-magnon spectral functions x^j (Co) 
for 5 = 1, 2 are shown in FIG. 9. For all values 

of 5, a clear two-magnon peak appears. In the most 
relevant case 5 = 1/2, its maximum is at w max — 1.29. 
Expressed in terms of the exchange coupling J , this 
translates into a frequency 07 max — 2.44 J, somewhat 
smaller than the result w max ~ 2.76 J obtained from 
an interacting spin-wave theory 23 . As will be shown 
below, the experimental spectra are consistent with our 
present result. Indeed the observed peak in the Raman 
spectra are close to 2.44 J with values of the exchange 
coupling which agree very well with those determined 
independently from neutron scattering. Taking into 
account the different frequency scales for different values 
of 5, the peak is shifted towards lower frequencies for 
smaller spin. This behavior has previously also been 
observed by Canali and Girvin 38 . A quite significant 
difference between their results and ours, however, is 
the pronounced asymmetry of the two-magnon peak. As 
will be shown in the following section, this asymmetry 
is essential for achieving a quantitative agreement with 
experiment. 












FIG. 9. (Color online) 2-magnon spectral function X 2 .Ma.gn 0 n 
for S — 1/2, S — 1 and 5 = 2. 


In order to see to which extent the value of the Higgs 
mass mo influences our results for the two-magnon re¬ 
sponse, we have varied this parameter in the range 0.1 < 
too < 0.9. Surprisingly, we find no detectable change in 
the two-magnon line shape. The physical reason for the 
fact that the two-magnon spectral function is essentially 
independent of the Higgs mass despite the fact that the 
magnon interaction is mediated by the amplitude mode, 
is the following: The inverse Higgs mass is the range 
of the interaction in position space while the interac¬ 
tion strength is controlled by the coupling g. As both 
magnons, involved in the interaction process, are created 
by a single photon, they are initially on neighboring lat¬ 
tice sites. For rh 0 < 1 the range of the Higgs-mediated 
interaction is larger than one lattice spacing, hence we do 
not expect a sensitive dependence of the 2-magnon peak 
on the Higgs mass. 

As a second step, we turn to the results for the two- 
Higgs susceptibility where a similar dependence on the 
spin quantum number S appears as for the two-magnon 
response. We therefore only show the results for S = 1/2 
which is relevant for the comparison with experiments 
below. In FIG. 10 the imaginary time data for the two- 
Higgs susceptibility is shown for different Higgs masses. 
First of all one notices, that the absolute values of the 
two-Higgs susceptibility are smaller than the ones for 
the two-magnon susceptibility. As a result, the two- 
magnon peak will always dominate the Raman spectrum. 
From the leading order analysis, we again expect a Cj~ 2 - 
behavior for large frequencies. An Ansatz of the form 
used in Eq. (33) again reproduces the numerical data 
rather well, with a variance 1.8 x 10 . The resulting 

conribution x^(cD) to the Raman spectrum in real fre¬ 
quencies is shown in FIG. 11. As expected, the spectra 



FIG. 10. (Color online) Numerical results for the 2-Higgs 
susceptibility as a function of imaginary frequency for mo = 
0.9, mo = 0.61 and mo = 0.25. The solid curves are fits of 
eq. (33). 



FIG. 11. (Color online) 2-Higgs spectral function Xh for S = 
1/2 and three different values of the Higgs-mass, mo = 0.9, 
mo = 0.61 and mo = 0.25. 


are shifted towards higher frequencies with increasing val¬ 
ues of the bare Higgs mass fho- In particular, the spec¬ 
tral weight is rather small below 2mo, reminiscent of the 
threshold behavior in leading order. 

An exact result which determines the relative weight 
with which magnons or the Higgs mode appear in the 
Raman spectrum can be obtained by considering the ra¬ 
tio of the first moments of the two-magnon and two- 
Higgs susceptibilities Xm and Xh- These moments can 
be expressed in terms of the exact progators G WT (fc) and 
G rJ(7 (k) in a way which is analogous to the standard f-sum 
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FIG. 12. (Color online) Ratio of the first moments, M 4 
and of the two-Higgs and two-magnon susceptibility, 

respectively versus the Higgs mass mo = mo/A in units of 
the momentum cut-off at leading order. 


rule for the dynamic structure factor. In explicit form, 
the sum rules are 



FIG. 13. (Color online) Total Ranran-spectrum, Xn amfm (k-0 = 
Xm(w) + Xh (w), for three different values of the Higgs-mass 
rho and S = 1/2. mo = 0.25 is an example for a light Higgs, 
mo = 0.61 for a intermediate Higgs and rho = 0.9 for a heavy 
Higgs. 


XfU) — / ™ 

~ / 27T 


did 


a4 1} = 


did 

2n 


w.XmM = 2 g J 7 (k) 2 G %n (k) 

(34) 

w.\h(w) = 2g J 7 (k) 2 G aa {k) 

(35) 


Using the propagators to zeroth order, the ratio be¬ 
tween the Higgs and the magnon part of the Raman re¬ 
sponse is controlled by the value of the Higgs mass mo- 
As shown in FIG. 12, the weight of the spectrum com¬ 
ing from the amplitude mode decreases continuously with 
increasing mass as expected on intuitive grounds. 

Finally, we consider the total Raman spectral function 
XRaman = Xm + Xh without the interference term, which 
is negligible in the regime g -C g c studied here. We again 
concentrate on S = 1/2 and show the Raman spectra 
for three different values of toq in FIG. 13. These val¬ 
ues correspond to a light, an intermediate and a heavy 
Higgs mass. For both a light and a heavy Higgs with 
m 0 < 0.35 or rho > 0.8, the full spectrum is dominated 
by the response associated with magnons. The contri¬ 
bution from x" only appears as an additional spectral 
weight at frequencies above the two magnon peak, while 
no distinct peak is associated with an excitation of the 
amplitude mode. It is only in the intermediate regime 
0.5 < rho 0.7 where a shoulder appears above the 
two-magnon peak as a remnant of the sharp onset of a 
Higgs-nrode in the leading order calculation of Eq. (20). 
As will be discussed below, such a feature is present in 
the experimental Raman spectrum of La 2 Cu 04 . 


E. Comparison with experiment 

In the following, we compare our theoret¬ 
ical results for the Raman spectral function 
X Raman with experimental data for the undoped 
cuprate materials YBa 2 Cu306 (Y-123) 13 and 
La 2 Cu 04 (LCO) 14 . Regarding the dominant two 
magnon response, the relevant energy scale is set by 
the characteristic frequency c s A of the underlying linear 
0(3) er-model, which is related to the exchange coupling 
by HJ = 7t _1 c s A 39 For undoped curates, the appro¬ 
priate values for J are known from inelastic neutron 
scattering. The agreement between the optimal values 
for J obtained from our fit to the Raman data and the 
n-scattering results is therefore a crucial consistency test 
for both the model and our approximations involved 
in calculating the Raman spectrum. Since neither 
experiments nor theory can reliably determine the 
Raman intensity in physical units of counts per second 
and per frequency interval, the results are plotted in 
arbitrary units. Correspondingly, the experimental data 
are fitted with a function of the form A\" (w/(c s A)), 
where an overall prefactor A remains undetermined. It 
is important to note however, that this arbitrary overall 
factor does not affect the relative weight to the Raman 
spectrum which comes from the magnons or the Higgs 
mode. Apart from the dimensionless Higgs mass rho this 
leaves the frequency scale c s A = 7 rhJ as the only physical 
parameter which the spectra depend on. An important 
and unexpected result of the comparison between theory 
and experiment below is that without the Higgs mode 
contribution the slow decay of the spectrum above 
the two magnon peak cannot be described properly. 
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The amplitude mode of 2D quantum antiferromagnets 
is therefore directly visible in the Raman spectrum 
despite the fact that generically it does not show up as a 
separate peak. It is only in the particular case of LCO, 
where the Higgs contribution gives rise to an additional 
shoulder in the tail above the two magnon peak which is 
likely due to an intermediate mass Higgs. 

The experimental data for the bilayer compound Y- 
123 are shown in FIG. 14 together with theoretical fits of 
the full Raman susceptibility. Apparently, the quite non¬ 
trivial form of the dominant two-magnon peak is repro¬ 
duced quite well by our theory with an exchange coupling 
J = 126meV which is rather close to the value 120meV 
obtained from inelastic neutron scattering 40 . The opti¬ 
mal value of the dimensionless Higgs mass for Y-123 is 
Too = 0.25. This is significantly smaller than the value 
mo — 0.6 which is found below for LCO. A possible origin 
for the small value of the Higgs mass, which entails large 
fluctuations of the magnitude of the Neel order, may be 
that Y-123 is a bi-layer system with an antiferromag¬ 
netic interlayer exchange coupling J±/J ss O.l 40,41 . Such 
a system undergoes a quantum phase transition from a 
phase of two weakly coupled 2D antiferromagnetic lay¬ 
ers to a gapped singlet phase with increasing J± 42-44 . 
While Y-123 is clearly in the phase where the two layers 
possess antiferromagnetic order, the interlayer coupling 
is expected to increase quantum fluctuations in the Neel 
order compared to the single-layer compound LCO and 
hence leads to a smaller Higgs mass. The presence of an 
amplitude mode directly shows up in the Raman spec¬ 
trum via the slow decay above the two magnon peak. 
Indeed, the dashed line FIG. 14, which is an optimal fit 
to the spectrum without the Higgs, clearly misses a sub¬ 
stantial part of the spectral weight in the range between 
3000cm -1 and 5000cm -1 . A feature which is not poperly 
described by our model is the slightly non-monotonic de¬ 
cay in this frequency range. Such a nontrivial structure 
may in principle be caused by a triple resonance discussed 
by Morr and Chubukov 45 . This resonance appears in sit¬ 
uations where the incoming photon energy fuvi is above 
the charge transfer gap 2A. The pure spin-photon in¬ 
teraction of the Fleury-Loudon Hamiltonian then needs 
to be extended by taking into account the generation of 
intermediate particle-hole pairs. While the experiments 
indeed were in a range with hu>i > 2A ~ 2eV, they do 
not show a change of the additional feature with increas¬ 
ing values of the incident photon energy, as expected for 
a triple resonance 45 . An explanation of the slightly non¬ 
monotonic decay of the spectrum in terms of a triple 
resonance therefore seems unlikely. 

For the single-layer compound LCO the data together 
with the theoretical fit are shown in FIG. 15. Again, the 
two-magnon peak is reproduced very well with a value 
J = 149meV for the exchange coupling which is in ex¬ 
cellent agreement with the INS result IdSmeV 1 '. Con¬ 
trary to the case of Y-123 and our theoretical result, the 
Raman spectrum of LCO exhibits a slight increase with 



FIG. 14. (Color online) Experimental data for the Raman 
spectrum of Y-123 together with a fit of the theoretical model 
with Higgs mass mo = 0.25 (solid line) and without Xh 
( dashed line). 



FIG. 15. (Color online) Experimental data for the Raman 
spectrum of LCO together with a fit of the theoretical model 
with Higgs mass mo = 0.595 (solid line) and without Xh 
( dashed line) 


frequency in the spectral range above 6000cm -1 . This 
deviation is most likely due to the luminescence of de¬ 
fects in the LCO sample, an effect which is absent in the 
Y-123 spectra shown in FIG. 14 because of a higher sam¬ 
ple quality 46 . Concerning the pronounced shoulder above 
the two-magnon peak at about 4500cm -1 , it turns out 
that this additional feature is reproduced quite well by 
choosing an intermediate value in o = 0.595 of the Higgs 
mass. The overall very good agreement between theory 
and experiment suggests, that the additional feature seen 
in the Raman spectrum of LCO is indeed associated with 
an amplitude mode of the underlying Neel state. What 
remains to be understood, however, is the fact that the 
pronounced feature the B\ g mode near 4500cm -1 in LCO 
also appears in the experimental data for the A 2 g mode, 
where no two-magnon response is visible. By contrast, 
there is essentially no Raman signal at any frequency in 
the A 2 g mode of Y-123 14 . 
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III. CONCLUSION AND OPEN PROBLEMS 

In summary, we have shown that the simple linear 
0(3) er-model is able to explain quantitatively the major 
features of the Raman spectrum of the Neel state of 
undoped cuprates. In particular, we have found that the 
detailed form of the dominant two-magnon peak is a re¬ 
sult of magnon-magnon interactions which are mediated 
by the Higgs mode. The relevant energy scale is fixed 
by the exchange coupling of the underlying Heisenberg 
model, with the peak located at w ~ 2.44 J. Apart from 
mediating interactions between the magnons, the Higgs 
mode also shows up more directly in an enhancement 
of the spectral weight above the two magnon peak. It 
leads to a separate peak only for an intermediate value 
of the Higgs mass, a case which is apparently realized 
in LCO. The agreement between experiment and theory 
depends quite sensitively on the proper inclusion of the 
Higgs mode contribution. Raman scattering therefore 
appears to provide a direct signature for the existence of 
a Higgs mode in 2D quantum antiferromagnets. 

While it is remarkable that a description based on 
the linear 0(3) er-model which only captures the low 
energy physics of quantum antiferromagnets is able to 
account for the detailed form of the observed Raman 
spectrum in the full range of energies up to 6 J and be¬ 
yond, our present study of possible signatures of a Higgs 
mode in 2D quantum antiferromagnets is certainly only 
a first step into a more detailed analysis of this problem. 
Among the open questions in this context we mention 
two: 

a) from a quite general point of view, what are suitable 
and experimentally accessible correlation functions 
where the Higgs mode in 2D quantum antiferro¬ 
magnets can be seen in direct form and can also be 
studied into a regime where Neel order disappears 
into some more complex spin liquid state? 

b) regarding the specific model studied above, can one 
calculate the mass parameter mo from an underly¬ 
ing microscopic model and - moreover- go beyond 
the perturbative calculation of the Raman spec¬ 
trum which is only reliable deep in the Neel phase? 

The complexity of the correlation functions which en¬ 
ter the Raman spectrum makes both problems quite chal¬ 
lenging. In particular, to uniquely identify the presence a 
Higgs mode from the experimental data so far, which in¬ 
clude the two cuprate examples discussed above and also 


the Iridate compound Sr 2 lr 04 25 , a better microscopic un¬ 
derstanding of the Higgs mass parameter is needed, both 
for single and double layer compounds, for which mo is 
apparently much smaller than in the single layer case. 
Regarding the issue of correlation functions where the 
Higgs mode shows up more directly than in a standard 
Raman experiment, an interesting option appears to be 
resonant inelastic X-ray scattering (RIXS) 47 . Of particu¬ 
lar interest here is indirect RIXS, which allows to measure 
a momentum dependent four-spin correlation associated 
with the operator 48 

Oq = J2 J k Sk- q • s„ k (36) 

k 

In the language of the NLcrM Eq. (36) corresponds to 
o q ~/d 2 x e* q x (d M n) 2 . In comparison to Eq. (9) one 
notices, that RIXS again couples to the square of field 
gradients, however in a symmetric form, without the an¬ 
tisymmetric tensor cr?-, which is responsible for the van¬ 
ishing of the leading diagrams associated with the Higgs 
mode. It is likely therefore that the Higgs mode in 2D 
quantum antiferromagnets is directly observable in a in¬ 
direct RIXS experiment, similar to the probe of a Higgs 
mode in the Bose-Hubbard model which is induced es¬ 
sentially by a modulation of the hopping amplitude 6 . 
By contrast, in direct RIXS one couples to a single spin 
operator (()G irect ) ^ ( ei x e/) • S q 49 . This is similar to de¬ 
termining the dynamic spin structure factor in inelastic 
neutron scattering and allows to measure the magnon- 
dispersion. 
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